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Quantum metrology is the state-of-the-art measurement technology. It uses quantum resources to enhance 
the sensitivity of phase estimation beyond what reachable within classical physics. While single parameter esti¬ 
mation theory has been widely investigated, much less is known about the simultaneous estimation of multiple 
phases, which finds key applications in imaging and sensing. In this manuscript we provide conditions of useful 
entanglement (among multimode particles, qudits) for multiphase estimation and adapt them to multiarm Mach- 
Zehnder interferometry. We discuss benchmark multimode Fock states containing useful qudit entanglement 
and overcoming the sensitivity of separable qudit states in three and four arm Mach-Zehnder-like interferome¬ 
ters - currently within the reach of integrated photonics technology. 


INTRODUCTION 

Quantum metrology exploits entanglement in the probe state to enhance the precision of parameter estimation beyond what 
is reachable with classical resources (see [1,2] for reviews). The role of entanglement in the estimation of a single parameter 
has been clarified [3-6] and investigated experimentally in Mach-Zehnder interferometers (MZI) [7]. However, much less 
is known about the role of entanglement in the joint estimation of multiple parameters. Yet, in many practical applications, 
multiple parameters are estimated simultaneously. This includes quantum imaging [8] as well as probing of biological samples. 
Interestingly, the theory of multiphase estimation does not follow trivially from what is known about the single parameter case 
[9]. Indeed, ultimate multiphase estimation bounds are not saturable in general [10], due to the non-commutativity of the 
operators generating the phase shift transformations [11, 12]. First insights on this scenario have been recently reported [13-18]. 

A natural platform for multiparameter quantum metrology is provided by multiport interferometry, generalizing conventional 
two-mode interferometry. Recent progresses in the realization of multiport devices have been achieved by exploiting integrated 
photonics [19-27]. Three- and four-port beam-splitters (tritters and quarters) have been produced with integrated optics [27- 
30]. This paves the way toward the realization of multiarm interferometers created by two tritters (quarters) in succession [31]. 
Quantum-enhanced single parameter estimation in integrated interferometers has been theoretically predicted [14]. 

In this manuscript we provide conditions of useful entanglement (among multimode particles, qudits) for the estimation of 
multiple phases and adapt them to the case of multiarm MZI. We analyze the simultaneous estimation of multiple phase shifts in 
an experimentally relevant framework, with multiphoton Fock states as probe and photon counting measurement. Our analysis 
generalizes the case of twin-Fock MZI which has attracted large experimental [7, 32-34] and theoretical [35-37] interest for 
quantum-enhanced single phase estimation. From the analysis of the Fisher information and employing an adaptive multiphase 
estimation, we predict a multiparameter estimation sensitivity beyond the limit achievable with separable qudit probe states. 


RESULTS 

Multiparameter estimation 

We consider here the estimation of a n-dimensional vector parameter A = (Ai, • • • , A n ) [9]. In our benchmark, every 
parameter corresponds to a phase to be estimated in a multiarm interferometer. A general approach (see Fig. 1 a) consists in 
preparing a probe state go, applying a A-dependent unitary transformation U\ and performing independent measurements on 
v identical copies of the output state g\ = U\QoU\. The measurement is described by a positive-operator valued measure 
(POVM), i.e. a set {Tl x } of positive operators satisfying J2 X IL X = 1, P(x\X) = Tr^nJ being the probability of the 
detection event x. Finally, the sequence x m (xi, • • • , x v ) of v measurement results is mapped into a vector parameter A (as) = 
(Ai (as), • • • , A n (x)), representing our estimate of A. A figure of merit of multiparameter estimation is the covariance matrix 

C»,j = Y P(X I A) [A* - Ai(*)] [Aj - Aj-(a:)], (1) 
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where P(x |A) = ]Xim P( x i |A) and A = (Ai, • • • , A n ) is the mean value of the estimator vector. For locally unbiased 
estimators (i.e. dAi/dXj = Sij) the covariance matrix is bounded, via the Cramer-Rao theorem [9], as 

c > F _1 /z/ (2) 


(in the sense of matrix inequality), where 


^ 1 dp(x | A) dp(x | A) 

1,3 “ P (a;| A) <9 A i d\j 


(3) 


is the Fisher information matrix (FIM). Notice that Eq. (2) can be derived only when the FIM is invertible. The equality sign in 
Eq. (2) is saturated, asymptotically in v, by the maximum likelihood estimator. In this Letter we quantify the phase sensitivity 
by the variance of each estimator, (SXj) 2 = Cjj. We have 


( a ,) 2 > 1 m > 


z/F 


3,J 


(4) 


where the first inequality is due to (2) and the second follows from a Cauchy-Schwarz inequality, see [38]. Since l/(z/F jj) is 
the Cramer-Rao bound for single parameter estimation, inequality (4) tells us that sensitivity in the estimation of A j is optimized 
when fixing the other parameters to known values. We will also consider 


!>,)*> ap > if 

j= 1 3 = 1 


(5) 


The right-hand side inequality in Eqs. (4) and (5) is saturated if and only if the FIM is diagonal, F ij = FiSij. Furthermore the 
FIM is bounded by the quantum Fisher information matrix (QFIM): F < Fq, where 

[Fq]»j = Tr [gxLiLj + Q\LjLi]/ 2, (6) 

and Lj is the symmetric logarithmic derivative of Q\ with respect to parameter A j, defined as djQ\ = ( LjQ\ + g\Lj)/2 [9]. 
In the single parameter case, the QFIM reduces to a single scalar quantity and it is always possible to find a POVM for which 
F = F q holds [39] and SX = 1 /Fq holds. In contrast, in the general multiparameter case, it is not possible to achieve this 
bound [10-12]. 


Sensitivity bounds for separable states 


Here we set sensitivity bounds for multiparameter estimation when the probe state is separable. A state go of N qudits is said 
to be qudit-separable if it can be written as g sep = AAkPkOk^ ® ® where (j = 1, • • • , TV) is a single qudit state. 

A state which is not qudit-separable is qudit-entangled. The notion of qudit (a particle in n > 2 modes) generalizes the concept 
of qubit (a two-mode particle) and is relevant in multimode interferometry [2]. We will set the conditions of qudit entanglement 
when the generator of each phase shift, Gj = is local in the qudit, i.e. it can be written as Gj = Y^iLi 9j^ where 

3 (i) 

acts on the ith qudit. For simplicity, we will take the same operator g - = cjj for each particle. For separable probe states the 

inequality 

F jj A(^j 5max ^jfjinin) (2) 


holds for all possible POVMs [38], where gj ^ max and g^min are the maximum and minimum eigenvalue of gj, respectively. 
Inequality (7) gives a bound on the sensitivity reachable with separable states for the estimation of the single parameter A j, when 
all other parameters are set to a known value. Inequality (7) can be always saturated by optimal measurements [38]. For the 
estimation of the single parameter, the violation of Eq. (7) is a necessary and sufficient condition of useful qudit entanglement [2, 
4]. Regarding the simultaneous estimation of multiple parameters, we can use Eq. (7) and the chain of inequalities (4) to obtain 


> 


N(9j ,max 


( 8 ) 


Inequality (8) is a bound of sensitivity in the estimation of the single parameter A j with separable states, when all the parameters 
are unknown. Summing Eq. (8) over all parameters, we obtain 




N •_! (9j ,max $ 7 ,min) 


2 ' 


(9) 
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According to Eqs. (8) and (9), for separable states such that the FIM is invertible, we recover - at best - the shot noise scaling of 
phase sensitivity, SXj oc TV -1 / 2 , which characterizes single parameter estimation. The optimal prefactor, which is equal to one 
in the qubit case, reflects the multimode nature of the problem. 


Multimode Mach-Zehnder interferometry 

In the following we discuss the estimation of a phase vector <p = (</>i,..., <p n ) in a multiarm Mach-Zehnder interferometer 
(MMZI) (see Figs. 1 b-c). The MMZI is built by cascading a d-mode balanced beam-splitter U^ - which is the natural 
extension of the standard 50-50 beam-splitter to more than two optical input-output modes [40] - a phase shift transformation 
Z7((/>) = being Ni the photon-number operator for the ith mode, and a second multiport beam-splitter U^ d \ 

Hence, this platform can be adopted as a benchmark to investigate simultaneous estimation of n = d — 1 optical phases. Indeed, 
it allows for a direct comparison between classical and quantum probe states, and can be adapted to represent a flexible platform 
for the analysis of multiparameter scenario by changing the unitary transformation of the input and output multiport splitters. 
The inequalities (8) and (9) thus read 


[F- J |,.i > k, and TrlF- 1 ] > (10) 

respectively. The recent experimental implementation of symmetric multiport beam-splitter [27, 30], by adopting integrated 
platforms, paves the way toward the future realization of optical MMZI. For d = 3 modes, U^ (tritter) has diagonal elements 
= 3 -1 / 2 and off-diagonal elements ( U^)i,j = 3 - 1 / 2 ^ 271-/3 w ith z ^ j. For d = 4 modes, (quarter) is (Z = 





FIG. 1. a. Conceptual scheme of multiparameter estimation, b. Three-mode MZI for two-parameter phase estimation built by two cascaded 
three-port beam-splitters. Phases (0i, 0 2 ) on modes (ki, k 2 ) are the parameters to be estimated, while (^ 1 , ^ 2 ) are two additional controlled 
phase-shifts c. Four-mode interferometer for two-parameter phase estimation built by two cascaded four-port beam-splitters. Phases (0i, (^ 2 ) 
on modes (ki, k 2 ) are the parameters to be estimated, while (0o, ^ 1 ,^ 2 ) are assumed known and controlled. 
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Three-mode interferometer, input 11,1,1 > 
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FIG. 2. a-c. Optimal phase sensitivity of the three-mode balanced MZI with 11,1,1) probe state and photon-number measurement. Contour 
plots of a. Tr[F -1 ], b. (F _1 )i,i, c. (F -1 ) 2 , 2 , as a function of 0i and 02. Tr[F _1 ] is minimized at the working points Q i and Q 2 (see main 
text), d-f. Optimal phase sensitivity of the four-mode balanced MZI with |1,1,1,1) probe state and photon-number measurement. Contour 
plots of d. Tr[F _1 ], e. (F -1 )i,i, f. (F _1 ) 2 , 2 , as a function of 0i and 02. These are shown for 0o — 0.001 to avoid undetermined points in 
the plot. The QCRB is achieved, for instance, at the point 0\ — [ 77 , tt\. Red areas indicate the violation of the separable bound defined in Eq. 
(10). Tables A and B report TrfF^ 1 ] and (F^ 1 )^ for different input states and their comparison with the separable bound (Sep). 


2 _1 and = —2~ 1 for i 7 ^ j. The phase vector is estimated from the measurement of the number of particles in each 

mode. 

As probe, we focus on multimode Fock states with a single photon in each mode. For the three-mode MZI, the probe is 
11,1,1), corresponding to the injection of a single photon in each input mode of the interferometer. The results of the calculation 
for F -1 are shown in Fig. 2 a-c. We observe that Tr[F -1 ] and the diagonal elements [F _ 1 ]i ? i and [F _ 1 ] 2,2 strongly depends 
on the phases 0i and 02 . Notably, the inequalities (10) are violated at certain optimal values of the parameters, signaling that 
the Fock state |1,1,1) contains useful qudit entanglement. Additionally, we observe characteristic features, (i) F 7 ^ Fq, in 
particular, the minimum value of Tr[F _1 ] is greater than the corresponding minimum value of the QFIM: minTr[F -1 ] = 
0.59 > Tr[Fg X ] =0.5 (see Fig. 2 a). This value is lower than the bound for separable states given by Eq. (10) and equal to 
Tr[F -1 ] > 0.667 (here N = 3 and n = 2 ). We also have min 0 1 ?( ^ 2 [F _ 1 ]i 5 i = 0.25 and min 0 1? ^ 2 [F _1 ] 2 ,2 = 0.25, which are 
smaller than the bound 0.33 for separable state (see Fig. 2 b-c). (ii) The FIM is not always invertible: at the phase values for which 
det F = 0 the bound ( 2 ) is not defined. Around these points (white regions in Figs. 2 a-c) [F _1 ]i i and/or [F _1 ] 2 ,2 diverges, (iii) 
The working points to obtain the minimum of the multiparameter bound do not lead to symmetric errors on the single parameters 
0i and 02 . More specifically, when Tr[F _1 ] = 0.59, the bounds for the error on the single parameters are different: 7 ^ 

( 502 1111 . This is obtained for instance for working point Qi = ( 0 i, 02 ) = (0.892,2.190), leading to ([F _ 1 ]i ? i, [F — 1 ] 2 , 2 ) — 
(0.282,0.309) and for working point Q 2 = ( 0 i, 02 ) = (2.190,0.892), leading to ([F _1 ]i 5 i, [F _ 1 ] 2 , 2 ) — (0.309,0.282), see 
Fig. 2a. 

In summary, with this choice of system and measurement it is not possible to saturate the quantum Cramer-Rao inequality 
simultaneously for the two parameters. Furthermore, according to point (iii) an adaptive estimation strategy (which we dis¬ 
cuss below) is necessary to obtain the minimum sensitivity on both parameters with symmetric errors, and thus saturate the 
multiparameter Cramer-Rao bound. 

We have repeated the above analysis for a four-mode interferometer (d = 4) with two unknown phases, 0 1 and 02 , and a 
known control phase 0o (see Fig 1 c). This configuration allows a comparison between three- and four-arm interferometers for 
the two parameter estimation. In the latter case the control phase 0 0 gives us an additional degree of freedom. We choose as 
input the Fock State |1,1,1,1). In Fig. 2 d-f the results of our calculations are reported for a fixed value of 0 0 , as well as the 
numerical analysis of det F. We observe that as in the previous cases the FIM depends on the value of the parameter to be 
estimated. Furthermore, also in the four-mode the achievtable sensitivity falls below the bound (10) for separable states: we 
have min0 1? 0 2 Tr[F _1 ] = 0.375, min0 1? 0 2 [F _1 ]i ? i = 0.1875 and min 0 1? ^ 2 [F _1 ] 2 ,2 = 0.1875 which are below the bounds 0.5 
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FIG. 3. Adaptive estimation of two phases, 0i and <p 2 with the three-mode interferometer injected by a |1,1,1). The adaptive protocol (see 
text) aims at reaching a phase uncertainty 50i ~ 5<p2 after v — 10000 independent measurements, a-b: Uncertainties Scpi/Scpm and 502/50 m 
obtained for different values of 0i and <p 2 (points) and normalized with respect to the expected value 50 m = 0.543/yT/ (see text). As an 
example, we report the results obtained for the specific cases (pi = n (c) and 02 = 7r (d). In these panels the blue line is the estimated value of 
(pi, the red line is the estimated 02. The inset shows the difference between the estimated value and the actual value of the phases, error bars 
are obtained by repeating p = 1000 times the numerical simulation of the protocol. 


and 0.25 given by Eq. (10) (N = 4 and n = 2, here), respectively. The most notable difference with respect to the previous case 
is that the QCRB is achieved, for instance in working point 0\ — [n, 7 r]. In addition, both diagonal terms are equivalent and 
only a two step protocol is needed (see discussion below). 

We have also compared the obtained results with the one achievable with other probe states. For instance, we consider a set 
of distinguishable particles |V^is) = ®q=i \q) (where | q) stands for a single photon on mode k g ), or an input coherent state | 0 ^) 
on input mode ki with a = \/3 for = 3 (a = 2 for <i = 4) and no phase reference. We obtain Tr[Fg X ] = 1 for both |V^ is ) 
and within the bound Tr[F -1 ] > 0.667 given by Eq. (10) for separable inputs. Similarly, Tr[Fg X ] = 0.75 for both |V4 is ) 
and |V4)> within the bound Tr[F -1 ] >0.5. Results are summarized in Table A and B. 


Adaptive phase estimation 

The above analysis has shown that the working regime in which the minimum uncertainty for the estimation of the two phases 
(pi and 02 with the three-mode interferometer does not give the same error on the two individual parameters. To overcome this 
limitation - and obtain approximatively a symmetric estimate of the two phases 0i and 02 - we exploited an adaptive algorithm. 
The protocol requires v independent measurements and the adoption of controlled phase shifts 0^ on modes k^, with i = 1,2, 
which have to be tuned during the protocol to perform the estimation in different working points (see in Fig. 1 b). It is divided 
in three-steps. In a first step, we set -01,2 = 0 and obtain a rough estimate of the phases (pi after a number of measurements 
much smaller than v. Then, in step 2 the tunable phases 0^ are adjusted so that (pi + 0* on arms 1 and 2 are set to be close to the 
working point Qi. In this step essentially half of the remaining resources are spent so as to obtain (0^^ + 0i) d= 50^^ and 
(02^ + 02) ± 502^ with an adequate estimator. Here (pf^ L \ 50represent respectively the estimation and the uncertainty 
of (pi around working point Q\. In step 3 the same procedure is repeated for working point Q 2 . Finally the tunable phases \h 1.2 
are subtracted so to recover 0 i , 2 ± 501 , 2 - 

The results are shown in Fig. 3 a-d where v independent measurements are performed. Half of the measurements (z/i = v/2) 
are performed in point Q 1 , where 50i ~ 0.531/y^ and 502 — 0.556/^/z^T, while the other half (z/ 2 = v/2) are performed in 
point Q 2 , where 50i 0.556 /a/F 2 and 502 — 0.531/^7^2. The expected error on a single phase 50^ after the two steps is then 

obtained as an appropriate combination of the values on the points Qi . More specifically, as the Fisher information is additive, 
the overall FIM reads F = zq_Fi + ^F 2 , where F^ is the FIM in working points Qi. We observe that the protocol permits to 
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achieve the bound of the working point, which for v\ = z/ 2 is £0i = £02 > 0.543/v^* Note that the bound is lower than the 
corresponding limit (10) for separable states £0* > 0.577 j\fv- 

The adaptive scheme for the four-mode interferometer is slightly different: in this case there is one working point, the point 
Oi, see Fig. 2, so we we apply a two-step protocol. In the first step, we obtain a rough estimate of the parameters with an initial 
error £. Then, in the second step we apply two supplementary phases Vh and to translate the working point of the protocol 
to the neighbourhood of 0 \. It should be noticed that a convergent estimation protocol in the second step requires to set 0o 
such that the quantity Tr[F -1 ] has no singularities. Note that the more 0o deviates from 0o = 0, the larger is the regular region 
around 0\ [38]. The price to pay is a slightly increasing the error in the estimation process. The results of the protocol for the 
four-mode case with 0o = 0.01 are then shown in Fig. 4 a-b. Similarly to the three-mode case, we observe that the protocol 
permits to achieve the bound of the working point, which is £0 i = £0 2 > 0.437/y^ = for 0 o = 0.01 (plane in Fig. 4), 
while the quantum Cramer-Rao bound reads £0^ > 0.433/y^- This shows that achieving a convergent numerical protocol leads 
to a slight decrease in phase sensitivity due to singular points in the neighbourhood of the working regions. 



FIG. 4. Adaptive estimation of two phases, 0i and 02 with the four-mode interferometer injected by a 11,1,1,1), for 0o = 0.01 and v — 10000 
independent measurements, a-b: Uncertainties and £02/£0™ obtained for different values of 0i and 02 (points) and normalized 

with respect to the achievable bound = 0.437/ yfu. The horizontal red line in the legend corresponds to the quantum Cramer-Rao bound 
for the single-parameter. 


CONCLUSIONS 

In this manuscript we have developed the general theory of quantum enhanced multi-phase estimation. In particular, we 
provide the conditions of useful entanglement among multimode particles (qudits) for the simultaneous estimation of multiple 
phases below the ultimate sensitivity limit achievable with separable states. In a realistic experimental scenario, using multi- 
mode Mach-Zehnder interferometers and photo-counting measurements, Fock state probes can be exploited for multiphase 
estimation with quantum-enhancement phase sensitivity. With respect to the estimation of a single phase, where Fock states are 
known to be a useful resource, our analysis evidences a much richer case: most notably, the phase sensitivity strongly depends 
on the phase value (the Cramer-Rao bound being not always definite) and on the interferometer configurations such as the three- 
and four-mode interferometers. Finally, we discuss and numerically simulate an adaptive estimation protocol which permits 
to achieve the expected bounds. The adaptive strategy becomes crucial in the multiparameter scenario since the simultaneous 
saturation of the ultimate limits for all parameters is in general not guaranteed. 

During the completion of this manuscript, a first implementation of a tritter-based interferometer for single-phase estimation 
has been reported [41]. 
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I. BOUNDS ON THE DIAGONAL ELEMENTS OF THE 
FISHER INFORMATION MATRIX 


II. UPPER BOUND TO 


Here we detail the demonstration of the inequality [1] 

- ¥~z (S1) 

To prove Eq. (SI) we recall that the square root of the positive 
definite Fisher matrix is given by the matrix with the same (or¬ 
thonormal) eigenvectors as F and the square root of its eigen¬ 
values. We indicate as fa > 0 and (with vjVj = 5ij) the 
eigenvalues and eigenvectors of F, respectively (Fi^ = fiVi). 
Notice that F is real and symmetric and thus diagonalize. In 
addition, F is positive semidefinite (therefore fi > 0) and as¬ 
suming that F is invertible we have fi ^ 0. We write 

F = J2fiVivJ, 

Vf = J 

(VF )- 1 = Vf= Y,- 7 t ViV i ’ 

i v/i 

and therefore 

Vf(\/f)“ 1 = VfVf-i = i, 

with 1 the identity matrix. Using the Cauchy-Schwarz in¬ 
equality 

1 = (^[v'FkfelVF 3 !]^) 2 

k 

k k 


we obtain 


1 < F iti [F-\i Mi, 


The following inequality holds 

Fi,i<Fg>[p(\j\, (S2) 

where F Q^[p(&)] — [ F o\h i is the diagonal element of the 
quantum Fisher information matrix, 

F^[m]=Tr[p( A)4 2 ] (S3) 


and Li is the symmetric logarithmic derivative, satisfying 


dp{ A) _ Ljp(\) + p(\)Li 

d\i 2 


The demonstration of the inequality (S2) follows from Ref. [2] 
(see also the review [3]). Using Eqs. (7) and (S4), we have 


F- • 

- 1 - 


ST' 1 f 0p(x|A) \ 

Ap(x|A)V d\, ) 


3?(Tr[p(A)r,n,]) 2 
V IV[p(A)n,)] 


where $l(y) is the real part of y. We then use the following 
chain of inequalities valid for all values of x: 

^{Tr[p(\)Lifl x ]) 2 < |Tr[p(A)L i fl,]| 2 

< Tr[p(A)n a; ]Tr[n ;t: Ljp(A)l/,], 

the first inequality due to ?ft{y) 2 < \y\ 2 and the second is due 
to Caushy-Schwarz. We thus obtain (for all values of x) 


Vt(TY[p(\)Lifl x }) 2 

Tr[p(A)fl,)] 


< Tr il x Lip(\)Li . 


(S5) 


Summing over x , we thus recover Eq. (S2) with f< q [pM] 
given in Eq. (S3). The equality sign can be saturated by taking 
a projective measurement n x on the eigenstates of Li. 

To find an explicit expression for Fq* [p(A)], let us write 
p( A) in diagonal form p( A) = |with pk > 0, 

{| k)} is a complete basis \k)(k\ = 1) and ^2 k Pk = 1- 
We have 


and recover Eq. (SI). The equality sign is saturated if and only 
if F = cl, where c is a positive real number. 


i$>[p(A)]=2 £ 

k,k' 

Pk+Pk'A o 


fa -p/c') 2 

Pk + 


(Ik\Hi\k ')\ 2 , 


(dim 
V d\i 


U~\ A)- 
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where 


Hi mi 


(S6) 



2 


In particular, 

F$*[pW] <4(AG,) 2 . (S7) 

The quantity [Fq\ i:i has the physical meaning of a single¬ 
parameter quantum Fisher information with respect to trans¬ 
formations (7(A) where all the parameters are fixed except A i 
(in other words = 0 for i ^ j). 

III. RELATION TO QUDIT ENTANGLEMENT 

The relation between [ Fq and qudit entanglement is a 
generalisation of the criteria of useful entanglement discussed 
in Ref. [4] (see also [5] and the review [3]). We recall that 
a state go of N qudits is said to be qudit-separable if can be 
written as 

&ep = X Pfe IV’mXV’mI ® ® ■■■\i>k,N){i>k,N\’ 

k 

where \fk,n) (n = 1, • • • , TV) is a single qudit state. We as¬ 
sume here that the Hamiltonian G{ is local in the qudits (i.e. 

it can be written as G{ = Yln=i 9^* where g acts on 

( n ) 

the nth qudit), and, for simplicity, we take g\ } = cp for all 
n = 1, • • • , TV. Under these conditions we find 

Fq* < A%i,max - 5i,min) 2 , (S8) 

where g^ max and g^ i n are the maximum and minimum 
eigenvalues of fa. Taking into account Eq. (S2), we obtain 
that the inequality 

Fi,i ^ TV(^,max Pi,min) (S9) 

holds for all separable qudit states and all possible POVMs. A 
violation of this inequality signals qudit entanglement. 

The demonstration of Eq. (S8) uses general properties of 
the single-parameter quantum Fisher information [3]): 

• the convexity of Fq* : 

Fq* [ft,ep] < X PkF, q* [IV’l.fc) ® ••• 0 IV’JV.fc)]; 
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• the additivity of Fq* • 

N 

F$* [\i>i,k) 0... 0 IVxfc)] = X f q [l^n,fc>]; 

n= 1 


• the bound 

F ( J ) [l^ )fc )]< 4 (Aft)X ,0 


• the inequality 

^(^Pi) \^ n — {pi ,max — Pi, min) 5 

which holds for every n, i.e. for every single-qudit state. 

Putting all these results together, we arrive at Eq. (S8). 

IV. THE ESTIMATOR FOR THE ADAPTIVE PROTOCOL 

The chosen estimator is essentially a Maximum likelihood 
estimator, with likelihood function £(</)): 

^max 

£(</>) = log P (a) (</>) H p(k\<f>) nk . (S10) 

k=l 

Here, p( a \f>) represents the knowledge on the parameters, 
p(k\cf)) is the conditional probability of outcome fc, and rik 
is the number of occurrence of outcome k. At each step, the 
distribution p^ (</>) is a Gaussian distribution: 

2 (4 >i — p 

P {a *(<t>) = ll , 1 e~ < 2 X )2 ) , (Sll) 

being and the estimated value and the error on the 
parameter i obtained at step a — 1 respectively. 
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SUPPLEMENTARY FIG. SI. Phase sensitivity of the three-mode balanced MZI with |1,1,1) probe state and photon-number measurement, 
a, Contour plot of Tr[F _1 ] as a function of fa and fa. b-c, Regions of Tr[F _1 ] around the values of the phases which minimize Tr[F -1 ]. 
d, working point Q i and c, working point Q 2 . d, Contour plot of (F _1 )i,i as a function of 0i and fa. e-f, Regions of (F _1 )i,i around the 
values of the phases which minimize Tr[F -1 ]. e, working point Q 1 and f, working point Q 2 . g, Contour plot of (F -1 ) 2,2 as a function of 0 1 
and fa. h-i, Regions of (F _1 ) 2,2 around the values of the phases which minimize Tr[F -1 ]. h, working point Q 1 and i, working point Q 2 . 



SUPPLEMENTARY FIG. S2. Analysis of stability for the four-mode 
interferometer with |1,1,1,1) probe state and photon-number mea¬ 
surement. To perform a two-step adaptive protocol for the estima¬ 
tion of two unknown phases with the highest precision, we need 
that in a neighbourhood of radius S centered in the working point 
0 1, the quantity Tr[F _1 ] has no singularities, a-b, Contour plot of 
Tr[F _1 ] around the working point 0\ — [tt, 7t] for a fa — 0.01 and 
b fa — 0.1. Green circles: regions with 6 = 0.01 around 0\. Red 
circles: regions with S = 0.1 around 0\. We note that increasing fa 
the singularities move away from the neighbourhood, increasing the 
stability of the estimation protocol. 



























































